On the possibility of superconductivity in PrBa2Cu307. 
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Recent reports about observations of superconductivity in PrBa2Cu307 raise a number of ques- 
tions: (i) of various theories striving to explain the Tc suppression in Pra;Yi_2,Ba2Cu307, are there 
any compatible with possible superconductivity in stoichiometric PrBa2Cu307? (ii) if this super- 
conductivity is not an experimental artifact, are the superconducting carriers (holes) of the same 
character as in the other high- Tc cuprates, or do they represent another electronic subsystem? (iii) 
is the underlying mechanism the same as in other high- Tc superconductors? I present an answer to 
the first two questions, while leaving the last one open. 



One of the most exciting cases of superconduc- 
tivity suppression in high- Tc cuprates is that of 
i?i!^i_a;Pra;Ba2Cu307, where RE stands for a rare earth 
(see Refs. for reviews). Even more exciting are re- 
cent indications that conductivity and superconductivity 
can be restored in pure stoichiometric PrBa2Cu307 
This is such an unexpected results ||] that it is still not 
generally accepted and further experimental confirmation 
is required. Nevertheless, this fact was reported by two 
independent groups, and it is time now to understand the 
theoretical consequences of this finding. The most impor- 
tant message, if this finding is true, is that at a: = 1, and, 
presumably, at intermediate a;'s, there are free carriers 
in i?£'i_2;Pra;Ba2Cu307, and the suppression of metal- 
lic conductivity at sufficiently large x must be due to 
localization of those carriers. This statement effectively 
eliminates the possibility of hole depletion due to hole 
transfer into occupied states ( "four-valent Pr model"). 
It furthermore becomes highly unlikely that any kind of 
magnetic pair breaking is in effect, because (1) normal 
state conductivity drops sharply with doping, indicat- 
ing the change of character, if not the number, of carri- 
ers, and (2) superconductivity is, supposedly, restored at 
X = 1. It seems that we can then consider only the mod- 
els which associate the (super) conductivity suppression 
with a transfer of holes to an itinerant, but different from 
the undoped YBCO, state, which should furthermore be 
prone to localization. At first glance the only theoretical 
model that satisfy this criterion is the itinerant model 
of Liechtenstein and Mazin [^|j7j . I will show below that 
contrary to the claim in the Fehrenbacher and Rice pa- 
per (FR), their model is also compatible with metallic 
and possibly superconducting behavior in PrBa2Cu307. 
1^ In fact, it turns out that the difference between the 
"local" FR model and the "itinerant" LM model is much 
smaller than it was thought to be; if handled correctly, 
the FR model never renders localized states. This new 
understanding means that the physics of superconduc- 
tivity suppression and its possible recovery is essentially 
the same in both model. There is still a quantitative 
difference between the two, which is hard to access ex- 
perimentally, but which is now of limited importance. A 
generic model that can be nicknamed FRLM explains the 
entire body of existing experimental results and does not 



seem to have any sensible alternative. An exciting fact is 
that this generic model not only provides a possibility for 
metallic and superconducting behavior of PrBa2Cu307, 
but that it also predicts the superconducting holes in it 
to be of entirely different physical nature than the carrier 
in familiar high- Tc superconductors. 

Let me start with a brief reminder of the essence of 
the FR model. The crystal structure of the YBCO 
family cuprates is such that a rare earth ion and four 
nearest oxygens form a nearly perfect cube. Moreover, 
among 7 orbitals of /-symmetry there is one, xyz, which 
has eight equivalent lobes directed along eight directions 
[±1,±1,±1]. In the standard coordinate system where 
X, y correspond to the CuO bond directions the same or- 
bital is {x^ — y'^)z. Since this /-orbital extends directly 
towards neighboring oxygens, one expects a noticeable 
pfa hopping between RE and O. Thus the electronic 
structure of a i?i?Cu204 bilayer breaks into two weakly 
interacting subsystems: usual Cu-0 pda bands, of which 
two antibonding ones cross the Fermi level in YBCO, 
and RE-0 fda states. Oxygen p states directed along 
the Cu-0 bonds ("per" orbitals) participate in the former 
and those perpendicular to the bonds ("p7r" orbitals) in 
the latter. If we start with a cluster of one Pr and eight 
surrounding oxygens and considered formation of an an- 
tibonding state of the f{x^-y^)z Pr orbital and eight oxy- 
gen p orbitals pointing directly towards Pr, we find one 
bonding, one antibonding, and 6 non-bonding states. If 
the energy difference between the bare O p level and bare 
Pr / level is not too large, the energy of this antibonding 
state, pfcr*, may become higher than that of the pda* 
Cu-0 state and will pull some holes out of the latter. 
Whether or not this will happen depends on the p — f 
energy separation and the p ^ f hopping integral. Sup- 
pression of superconductivity in i?_Ei_a;Prj:Ba2Cu307 is 
thus ascribed to the hole transfer from the superconduct- 
ing pda* band into the pfa* state. An indispensable 
component of this model is localization of carriers pro- 
moted into the pfa* state. FR argued that the oxy- 
gen orbitals forming this state form the 45° angle with 
the CuO planes and thus the orbitals of the same oxy- 
gen pointing towards neighboring Pr ions are orthogonal 
to each other. They also neglected direct hopping be- 
tween the oxygen pn orbitals. In such an approximation 
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the effective bandwidth of the pfa* band is zero, and 
the holes there are localized by infinitesimally small dis- 
order. This was the original explanation of the lack of 
(super)conductivity in PrBa2Cu307. Furthermore, since 
this model renders noticeable presence of Pr states at the 
Fermi energy, one expects the Curie temperature for or- 
dering of the Pr moments to be much higher than for 
other _Ri?Ba2Cu307, which is indeed the case. A draw- 
back is that total localization in the FR model does not 
let the given PrOg cluster be influenced in any direct 
way by the rare earths filling other cells, in contradiction 
with the experiment: the Tc suppression rate even at low 
doping depends strongly on the host rare earth. 

The LM model Q differs from the FR model in es- 
sentially only one aspect: direct hopping between oxy- 
gen orbitals jl^ is taken into account. LM calculated 
this hopping as well as other relevant parameters of the 
electronic structure numerically using LDA-I-U method 
including Coulomb correlation in the rare earth /-shell. 
They found substantial 0-0 hopping so that independent 
of the value of the pfa hopping (and even the presence 
of the / orbital) the FR states were forming a dispersive 
band. It was originally thought ^ that an advantage 
of the LM model over the FR model was that a disper- 
sive band would hybridize with all rare earth ions in the 
crystal, and its position before doping with Pr would de- 
pend on the position of the /-level in the host rare earth. 
This naturally, and with reasonable quantitative agree- 
ment, explains the different rates of suppression with dif- 
ferent host RE H]. This model was however criticized 
[ pT[ because localization of carriers in such a dispersive 
band requires finite disorder and one expects stoichio- 
metric PrBa2Cu307 to be metallic. Another prediction 
which was seemingly different from that of the FR model 
was that the holes transferred to the pfcr* band con- 
centrate near the center of the Brillouin zone; that is, 
near k= (7r/a,7r/6). This prediction could be indirectly 
checked by measuring the ratio of the out-of- plane (pz) 
and in-plane (px.y) oxygen characters. At least qualita- 
tively, this prediction was confirmed by the experiment 

Interestingly, it was not noticed until very recently 

that the geometric argument of FR j|] was incor- 
rect: in fact, the angle that an 0-Pr bond forms with 
the CuOa planes is not 45°, but tan-i(l/V2) « 35°16', 
which means that even in the FR limit of no direct 0-0 
hopping, the pfa* states form a band whose dispersion 
is defined by the Pr-0 hopping amplitude tpf^. Below 

1 show how this band forms, using the nearest-neighbor 
tight-binding Hamiltonian. 

Let us begin with some notations: first, we neglect the 
(very small) z-dispersion. (This means that all orbitals 
we consider are antisymmetric with respect to z ^ —z 
reflection, like the fz{x^-y^) orbital.) Then the two plane 
problem is equivalent to a single plane. If we include 
all nearest neighbors, the following orbitals contribute 
to the FR band: (1) Pr z{x^ - y^), (2) 02 z, (3) 03 
z, (4) 02 y, (5) 03 x, (6) Cu xy, (7) Cu yz, and (8) 



Cu zx. Their in-plane 2D symmetries are, respectively, 
— y^, s, s, y, X, xy, y, and x, which simplifies the 
task of the tight-binding description of the band struc- 
ture. Let us now identify the largest hopping amplitudes 
between these orbitals. According to FR, this is pfa, 
which we shall denote tpf. It controls the following hop- 
pings: ti2 — tl3 = \[^^pf^ ^14 = ii5 = \J^^pf — V2ii2- 
This parameter defines the effect of the rare earth sub- 
stitution on the FR band. Another hopping, which is the 
strongest according to LM, is of pdi: type, denoted tp^. 
The hopping amplitudes ^28. ^37. ^461 and ^sg all are equal 
to tpd- This parameter defines the dispersion of the FR 
band in the absence of the / states; e.g., in YBa2Cu307, 
or in the spin-minority channel of PrBa2Cu307. Let us 
first consider these two hoppings separately. 

The FR model corresponds to an approximation tpd — 
0. Dispersion of the oxygen pir states is completely ne- 
glected; an isolated Pr impurity forms a localized anti- 
bonding state, shifted up with respect to the bare O p 
level by 



ek 



5e 



9Er,-E 



f 



(1) 



where 8 stands for the eight neighboring oxygens. It is 
assumed that tpf <^ Ep — Ef. In the opposite limit, when 
all rare earth sites are occupied by Pr, a narrow band is 
formed with the dispersion 

fk — -Bp = (5e — Secos2ip{cosakx + cosbky)/2, (2) 

where ip = arctan(l/-\/2) is the angle that the Pr-0 bond 
forms with the xy plane. Had this angle been 45°, as 
assumed by FR, the band would be dispersionless and 
thus fully localized. In reality, it should acquire a finite 
bandwidth W = 5eco82ip — Se/3, even were without 
Pr. This is an example of dispersion due to nonorthog- 
onality: the Hamiltonian written in terms of the oxygen 
orbitals pointing towards Pr is diagonal, but such a ba- 
sis is nonorthogonal and that results in dispersion. Note 
that the top of the band occurs at the (tt, tt) point and 
that is where the holes go from the pda superconducting 
band. Fig.|l| illustrates that indeed at this point the pfa 
interaction is antibonding along all bonds. 

Now we consider the case of finite tpd and no / states. 
For simplicity, we let the energy of the Cu d level be 
the same as the energy of the O p level. Then four O p 
orbitals and three Cu d orbitals form three antibonding 
bands (besides the bonding and nonbonding bands): 
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The top of the highest (third) band is again at (tTjTt), 
as iUustrated in Fig. showing again the antibonding 
interactions for all bonds. 

The case of pdn and pfcr interactions taken together 
cannot be solved analytically. Before reporting the nu- 
merical results, we make one additional observation: 
since both cases separately produce dispersive bands with 
the maximum at (7r,7r), one might expect this effect 
(band dispersion) to be enhanced when both interaction 
are included. One can easily see (Fig.p, that this is not 
true: the configuration of the O p orbitals, which is an- 
tibonding in the first case, is nonbonding in the second 
case, and vice versa. 

The main shortcoming of the original FR model was 
its inability to describe the different rate of Tc suppres- 
sion with the different rare earth hosts. The LM model 
with its dispersive pfa provides a natural explanation. 
However, we can see comparing Eqs.(||) and (H), that at 
low doping, that is near the (tt, tt) point, the shape of the 
FR band is very similar to that of the LM band. Note 
that the scale of the dispersion, i.e. the effective masses, 
may be different — the FR band should be heavier than 
the LM band and thus easier to localize. However, for 
the ideal stoichiometric PrBa2Cu307 at zero temperature 
both models give a metal and possibly a superconductor. 

Another indirect argument in favor of the LM model 
over the FR model was deduced from recent near-edge 
X-ray absorption experiments jl2|. ft was found that O 
Ptt orbitals form a relatively small angle with the Cu02 
planes: 20 to 25° at the doping level of 80% Pr. This is 
closer to the prediction of the LM model (15 to 18°) 
than to the original FM prediction (45°), and even to 
the corrected number of 36°. However, the fact that the 
correct FR model is nonorthogonal not only yields a fi- 
nite bandwidth in pure PrBa2Cu307, but also makes the 
average angle of the O orbitals dependent on doping 
at small doping — similar, but quantitatively different 
from the prediction of the LM model. 

Indeed, the orbitals around an isolated Pr impurity 
are tilted hy ip = arctan(l/'\/2); it is however obvious 
that if an oxygen atom has Pr ions on both sides, the 
highest state does not include 0(pz) character of this 
oxygen at all. In the low doping limit, for the Pr concen- 
tration X, the probability for an oxygen to have 1 Pr 
neighbor is vi = 2x{l — x), and to have 2 neighbors 
is V2 = x^ . Thus the average p^ character for the in- 
plane oxygen holes, seen in an experiment like Ref. | [T^ ] 
is Uz = vi SVC? (p, while the total number of the holes in 
the FR state \s n = vi + V2. The average tilting angle 
is thus sin^a ~ f^Ef- For x = 0.8 we find a « 19.5°, 
in excellent agreement with the experiment. Moreover, 
this number is the lower bound on a in the FR model, 
because at large x one cannot neglect dispersion of the 
FR band, which will force a to deviate from the formula 
above (at x = 1 the FR model should give the same num- 
ber as the LM model, which is about 20°), so that at 
X — 0.8 we find a ^ 20°. Interestingly, while both the 
FR and LM models predict a dependence of the angle 



on concentration, and both must give the same value at 
X — 1 [a does not depend on the effective mass at this 
point), they predict the opposite dependences: in the LM 
model a falls to zero when x — > 0, while in the FR model 
it increases up to a = « 36°. 

My conclusion is that after being corrected to take into 
account the right geometry of the Pr-0 bonds, the FR 
model provides better agreement with the experiment 
than the LM model. I emphasize that after such cor- 
rection the difference between the two models is not the 
difference between a band model and a localized model, 
but the difference between two band models, one where 
dispersion originates from the pfa Pr-0 hopping, and 
another where it appears mostly due to the pdi: Cu-0 
hopping. The fact that the former appears more success- 
ful (however, the final word will be said by an experiment 
accessing x dependence of the angle a), does not mean 
that the pdn hopping is negligibly small. As illustrated 
on Fig.lc, it merely means that the dispersion due to this 
hopping is weaker than that due to ip/o-. 

To summarize, current experimental situation in 
Pr3,i?£'i_2,Ba2Cu307 is such that the band version of 
the Fehrenbacher-Rice model presented here explains 
all existing experiments addressing superconducting and 
transport properties of this system, including the recent 
observation of superconductivity at full substitution. In 
fact, PrBa2Cu307 is a more novel superconductor than 
all other cuprate high materials known: it is the only 
one where superconducting carriers are not residing in 
the Cu(a;^ — y^)—0{pcr) bands, but are of entirely differ- 
ent character. 
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FIG. 1. Tight-binding orbitals at the point S=(7r,7r) pro- 
jected onto X — y plane. Upper paneh FR model, tpd = 0. 
The Pr ion is in the center. Middle panel: LM model for 
YBa2Cu307, no Pr z{x^ — y^) orbital, finite tpd- Cu2 ions are 
in the corners. Lower panel: Illustration of inability of the 
02 y and 03 x orbitals to make an antibonding combination 
simultaneously with the Cu2 xy orbital (in the corners) and 
Pr / orbital (in the center). 
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